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Abstract 



Isospin violating signals in the r" — > (37r)^^'r decay mode are discussed. 

^ . For the t~ — > Tr^vr^vr+z^T- decay mode, isospin violation arises from the vector 

QQ , current contribution in the t~ -^ ujtt"^^ decay with the subsequent isospin vi- 

C^ I olating Lo decay into tt+tt". We demonstrate that such effects may be observed 

f^ I in presently available data through the measurement of the interference effects 

t^^ ' of these vector current contributions with the dominating axial vector current, 

5^^ ■ i.e. through a measurement of the structure functions Wp, Wg, Wh and Wj. 

'"^ In the case of the t~ -^ ir^-K^TT^v-r decay mode, a vector current contribution 



i-C ' ence terms being too low for present statistics. 



is generated by r]iT^ mixing in the decay chain r -^ rjp Ur ^ vr'^vr'^Tr v-^. We 
Q I find that this effect is rather small, the magnitude of the associated interfer- 



I. INTRODUCTION 

Isospin rotations have been successfully used in r decays into an even number of final 
state pions to relate the vector current to the corresponding cross sections measured in 
electron positron collisions ^^- In the case of the two pion mode, the r decay rate has 
been measured with a relative error below one percent which is of the size of possible isospin 
violating effects. Isospin symmetry relations are also very useful to relate various decay 
amplitudes in SttUt-, KKirUj. and KtitiVt- final states 0-§|]. 

Isospin violation effects in the decay r~ -^ utt'i^t have been discussed in [|]. Such signals 
could be revealed by an analysis of the angular distribution in the ujtt' system. Another 
interesting isospin violating process is provided by the decay r^ -^ rjn^Ur- The different 
theoretical predictions for the branching ratio [[^Q are still one order of magnitude smaller 
than the actual experimental upper limit 0. 

In this article we will concentrate on possible isospin violating effects in the r -^ Stti/^ 
decay mode. Although the theoretical uncertainties in this decay mode are fairly large, 
observations of small isospin violating effects (below 1% to the rate) might be possible 
with presently available statistics. The sensitivity to such small effects is provided by an 
analysis of angular distributions. The relevant information is encoded in structure functions 
||TI1|,TT] which allow to reconstruct the form factors in the dominating axial current and 



in the small isospin violating vector current contributions. In particular the interference 
effects between the vector and the axial vector amplitudes, given by the structure functions 
Wp, Wgi Wh, Wi allow for such a measurement. Any nonvanishing contribution to these 
structure functions would be a clear signal of isospin violation in the three pion decay 
mode of the r. After specifying the isospin violating vector form factor, we will present 
numerical predictions for these structure functions including the full dependence on the 
resonance structure. We also analyze Dalitz distributions for the purely axial vector structure 
functions. Such distributions, in particular for the structure function Wd, are fairly sensitive 
to the details of the p sub-resonance implementation in the underlying models. 

A branching fraction of 0.6% in the r — ^ Sttz/,- mode due to isospin violation has been 
reported by the ARGUS collaboration |]T2[. Their analysis is based on a study of eight 
different contributions to the amplitude. Unfortunately the relevant interference terms with 
the isospin conserving part of the amplitude cannot be traced out unambiguously from that 
work. 

The paper is organized as follows: The general structure of the decay amplitude and 
the structure function formalism in the three meson decay mode is briefly summarized in 
Sec. H and a particular choice for the form factors in the axial vector current, the Kiihn- 
Santamaria model [Q], is specified in Sec. |T|. Isospin violating contributions to an additional 
vector current form factor will be discussed in Sec. |3 (Sec.0) for the r~ — > 7r~7r~7r"'"i/r 
(r~ -^ TT^TT^TT^Ur) dccay mode. The relevant hadronic matrix elements are determined in the 
vector meson dominance model. We obtain as a by-product from the decay r^ -^ rj-K^Ti^v^ 
a new parametrization for the transition of the p resonance into three pseudoscalar mesons 
which is also needed as an input e.g. for the r decay into KKtcUt final states. Finally, 
isospin violating signals induced by the vector current form factor are discussed in Sec. VI. 



II. THREE MESON DECAY MODES: 
FORM FACTORS AND STRUCTURE FUNCTIONS 

The matrix element Ai for the semi-leptonic r decay into three mesons hi,h2, h^ 

t{1,s) -^ frills') + hi{qi,mi) + /i2(g2,"^2) + h^{qz,m^) (1) 

can be expressed in the following form: 

M{r^v^ hh,h) = ^ C';) u{l', s'h,{l - 75)«(/, s) .F. (2) 

In Eq. (0) Gp denotes the Fermi-coupling constant and O^. is the Cabibbo angle. The hadronic 
current 

J'{qi,q2,q,) = {h{q,)h2{q2)hiqs)W{0) - A^^mO) (3) 

is characterized by four independent form factors ^1,^2,^3, F4 [|l^. These form factors are 
in general functions of si = (g2 + 53)^,52 = (?i + q-sY, S3 = (gi + ^2)^ and Q^, which is 
conveniently chosen as an additional variable. 

J^iqi, 92, qs) = Vr Fi + V^/ F2 + t Vi" Fs + V^ F4 (4) 

with 



vr 


= {qi-q3)uT^'' 


vr 


= {q2-q3)uT>^'' 


vr 


= e^"'^^gi„g2/3g3 7 


vr 


= q^ + qll + q'3 = 



(5) 
T'^'^ denotes the transverse projector 

Tf,u = g^u - —^ ■ (6) 

Fi and F2 determine the spin one component of the axial vector current induced amplitude, 
F4 the spin zero part which is given by the matrix element of the divergence of the axial 
vector current. The vector current induced amplitude is responsible for the form factor 
F3. All form factors may contribute in the general three meson case [^,0,^. G-parity 



conservation and PCAC in the three pion decay mode implies ^3 = ^4 = 0. However, 
isospin violation is expected to give a nonvanishing contribution to F3. Such contributions 
will be studied in the last three sections of this paper. 

The three meson decay in Eq. ([1|) is most easily analyzed in the hadronic rest frame 
qi+q2 + qs = 0- The orientation of the hadronic system is in general characterized by three 
Euler angles (a, P and 7) as introduced in |IO| , p^ . Of particular interest are the distributions 



of the normal to the Dalitz plane and the distributions around this normal. Performing the 
analysis in the hadronic rest frame has the advantage that the product of the hadronic 
tensor {H^'^ = 7^(7^^)^) and the leptonic tensor reduces to a sum L^^H^^ = J2x ^xWx- 
The leptonic factors Lx factorize the dependence on the Euler angles. For the definition of 



these angles and the exphcit dependence of the coefficients Lx on a, /3 and 7 see ref. [|T0[- The 
(in general 16) hadronic structure functions Wx correspond to 16 density matrix elements 
for a hadronic system in a spin one and spin zero state (nine of them originate from a pure 
spin one state and the remaining originate from a pure spin zero state or from interference 
terms between spin zero and spin one). These structure functions contain the dynamics of 
the three meson decay and depend only on the form factors Fi and on the hadronic invariants 
Q^ and the Dalitz plot variables Sj. The scalar contribution is expected to be small |T^ for 
all three meson final states and will be neglected in the subsequent discussion of this papeiQ. 
Instead of the 16 real structure functions which characterize the general hadronic tensor 
H'^'^ one thus deals only with nine functions Wx- These nine structure functions can be 
divided in four functions which arise only from the axial vector current {Wa,c,d,e), one from 
the vector current (Wb) and the remaining four from the interference of the axial vector 
and vector current {Wf^g,h,i)- The latter will be of particular importance in the subsequent 
discussion. 



The dependence of the structure functions on the form factors F^ reads |T^ 
Axial vector structure functions: 

Wa = {xl + xl) \F^\^ + {xl + xj) IF2P + 2(xiX2 - xj) Re (F^F;) 
Wc= ixl-xl)\Fi\'' + {xl-xl)\F2\^ + 2{xiX2 + xl)Re{FiF;) (7) 

Wd= 2 [0:1X3 |Fi|2 - X2X3 IF2P +X3(X2 - Xl) Re(FiF2*)] 
We = -2x3{xi + X2)lmiFiF;) 
Vector structure function: 

Wb= xllF^l^ (8) 

Axial vector-vector interference structure functions: 

Wf = 2x4 [xi Im (F1F3*) + X2 Im (F2F3*)] 
Wg = -2x4 [xi Re (F1F3*) + X2 Re (F2F3*)]] 

Wh = 2x3X4 [ Im (F1F3*) - Im (F2F3*)] (9) 

Wj = -2x3x4 [Re (F1F3*) - Re (F2F3*)] 



The variables Xi are defined by 



Xl = V^ = ql-ql 

X2 = V^ = q^ - ql 

X3 = Vr = ql = ~ql (10) 

X4 = V^ = jQ^Xsq^ 



^Using an ansatz for a scalar contribution in the SirVr decay mode as specified in |T^, U. Miiller 
constrained such a contribution in the branching ratio to be less than 0.84 % by analyzing the 
spin- zero-spin-one structure functions with 1994 OPAL data [16|. Note that a possible scalar 



contribution would not contribute to the vector-axial vector interference structure functions in 
Eq. (y) which are important for an observation of isospin violating effects. 



where qf (qf) denotes the x (y) component of the momentum of meson i in the hadronic 
rest frame. They can easily be expressed in terms of Si, S2 and S3 [0]. iy^(Q^,Sj) and 
WsiQ"^, Si) govern the rate and the distributions in the Dahtz plot through 

The remaining structure functions determine the angular distribution. All of them can be 
determined by a measurement of the (3 and 7 dependence even without reconstructing the 
r rest frame. 



III. AXIAL VECTOR CURRENT CONTRIBUTION TO r" -^ {3ir)-i^r 

T decays into three pions are dominated by the axial vector current which allows for sig- 
nificant simplifications: G-parity implies F3 = 0, Bose symmetry relates Fi and F2 through 
^2(^1, S2, Q^) = Fi{s2, Si, Q"^) and PCAC leads to F4 = 0. Note that the structure functions 
Wb,f,g,hj in Eqs. (||g) vanish for F3 = 0. 

The two like-sign pions in r" -^ ir^ ir^ n^ Ur and r^ -^ tt'^it^tt^i't- are labeled such that 
\p2 1 > |pi I and p3 refers to the unlike-sign pion. The normalization of the form factors Fi 
and F2 for the three pion decay mode is determined in the chiral limit^ |p!7| , 

Fi = F2 = t'^, U = 93 MeV (12) 

For large Q^, Si and S2 these form factors are modulated by resonances in the Svr and 2tt 
channel. Following the ansatz of Kiihn and Santamaria |l[], one has 

Fi{Q\s2) = v^ BWa,{(f)Tf^\s2) (13) 

The Breit-Wigner functions BWx(s) are parametrized including energy dependent widths, 

2 
777 

m^x{s) = — ""■ rv i \ ' BWx(0) = l (15) 

rrix- s -i^sVx{s) 

For the a^ we have in particular 

Va,{s) = ^ F,, -f^ ma, = 1.251 GeV , F,, = 0.475 GeV (16) 

where the function g{s) has been calculated in [|l| and is derived from the observation, that 
the axial vector resonance ai decays predominately into three pions. 



^We use the Condon-Shortley phase conventions. 
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FIG. 1. Invariant mass distributions x 



/Q^ = m{Tr~iT~Tr~^) (solid), x = y/si^ = m{TT~^iT~) 
(dashed) and x = y/s^ = m{iT~TT~) (dotted) of the structure functions Wa (a), Wc (b), 
sign{si — S2)Wd (c), sign{si — S2)We (d) in the t~ -^ (Stt)"?/,- decay mode. 

The superscript (2m) in the p form factor T^'^"^\s) denotes the subsequent decay into two 
pions. In the parametrization of T^^'"^(s) one allows for a contribution of the first excitation 
P', 



rp{2m)/^ 



l+Pp 



BWp(s)+/5pBWp,(s) 



(17) 



with the energy dependent width 

'' '^ s yrrij^ — Ami J 

and similarly for the p'. The parameters are given by 



(18) 



6 



/3p = -0.145 



rup = 0.773 GeV , Tp = 0.145 GeV 



rrip' 



1. 370 GeV, r^/ = 0.510 GeV. (19) 



which have been determined from e~^e~ —>■ tt+tt" in |jl|. Predictions for the (si — S2)- 
integrated structure functions wxiQ"^) = J dsids2Wx{Q'^,Si,S2) for X = A,C,D,E based 
on this model are in good agreement with data [|1^. The invariant Stc and 2tt mass 
distributions for the four integrated nonvanishing structure functions WA,Wc,sign{si — 
S2)WD,sign{si — S2)We in Fig. [l| reveal the importance of the ai (solid) and p (dashed) 
resonances. The y/s^ distribution (dotted line) is then fixed by phase space restrictions and 
the y/Q^ and y/s^ ^ distributions through s^ = Q'^ — Si — S2 + 3m^. The structure functions 
Wd and We are combined with an energy ordering sign{si — S2) to account for Bose symme- 
try. The y/s^ 2 distributions of Wa, Wc and sign{si — S2)Wf have a clear peak around the 
p resonance, whereas sign{si — S2)V^z)(v^i 2) has a surprisingly different behaviour in the 
Kiihn-Santamaria model. The distribution shows a relatively wide peak around ^/s^ ^ — 0-5 
GeV and only a much smaller additional peak around the p mass. In contrast, the ^/s^ 2 



distribution for sign{si — S2)Wd based on the model in |T^ has its maximum around the p 
mass and only a small additional peak around ^/s^ ^ — O-^- An experimental confirmation 
of the predictions for the y/s^ ^ distributions in the axial vector structure functions shown 
in Fig. m and in particular in sign{s\ — S2)Wy) would be a good test of the details in the p 
resonance structure in the Kiihn Santamaria model which we use for the two axial vector 
form factors F\ and F2. 

IV. VECTOR CURRENT CONTRIBUTION TO r -^ vr vr 7r+:^^ 

More detailed studies, such as testing the magnitude of amplitudes induced by F3 through 
isospin violation, are possible and will be discussed in the following. Since they affect the 
angular distributions through interference terms between the (small) contribution from F3 
with the large contributions from F\ and F2 (see Eq. (H)), they should be accessible in 
measurements of the structure functions Wf,g,h,i, already with the statistics of ongoing 
experiments. 

A small vector current contribution (~ F3 in Eq. (^)) to the r^ — > tt^tt^tt+i/^ mode 
is expected to arise from the r~ —>■ wn^v^ decay with the subsequent isospin violating u) 
decay into tt+tt". G-parity requires that the tuvr^ system is in a 1~ state and hence the 
T^ -^ ujTi^Ur decay can only proceed via a vector current. The hadronic matrix element is 
determined through [0,|3l (^^ another approach see ||2T[| ) 

(a.(gi,A)7r-(g2)|F^(0)|0) = te'^'^^'' s*M, X) qif^q2, Ft\Q') 

Fi,^^\Q^) = I IfZ^ Ti^")(Q2) (20) 

nip ^ 



where V^ is the vector part of the weak current. Note that we have fixed the sign of Fy {Q) 
from 7r° — > 77. L- is the coupling of the charged p^ to the gauge boson W^ and is related 
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FIG. 2. The invariant ujtt mass spectrum from r -^ ujtt v^ measured by ARGUS [|2^ (filled 
circles), by CLEO p^] (opened circles) and by ALEPH ||25[ (opened squares). The sohd line shows 
the fit result to Eqs. (|o|,||). 



to the p°7 coupling /p, gp^T, is measured in the decays uj 
respectively, 



7r°7 and lo 



TT TT TX 



22 



/p- = 


y2/p ~ 0.17 GeV^ 






y pujTv 


' 11.7 ±0.4 GeV"^ 
15.0 ±0.1 GeV"^ 


from T {uj - 
from T iuj - 


-.vrV) 



(21) 



The p-meson and its radial excitations are possible resonance candidates for the vector 
form factor T(^™^(s), where the superscript [Am) refers to the (anomalous) VMD decay 
chain p -^ uj-n -^ An. The admixture of the radial excitations in T(^™'''(s) is expected to 
differ from the corresponding p form factor To^^™-'(s) with a dominant two pion decay in 
Eq. (p!?!). Here we allow for an admixture of the p' and the p" via 



Tj^™)(s) 



1 



l + \ + K 



BWJs) ± ABW„/(s) ± KW^,n{s] 



(22) 



where we fix the parameters to the PDG |^2[ values, which yields 



mp = 0.773 GeV , 

TUp, = 1.465 GeV , Tp, ■- 

nipii = 1.70 GeV , Fp// 



0.145 GeV 
= 0.310 GeV 
= 0.235 GeV . 



(23) 



The parameters A and /t are obtained from a fit to the normahzed invariant mass spectrum 



of r -^ ujTT Uj. data P3|-p5[], see Fig. g, 

A = -0.054 ±0.012 

n = -0.036 ± 0.004 (24) 

xVd.o.f. = 53.3/31 . 

Note that the errors should be taken as an educated guess only, since we fit to the published 
data with the correlation matrices to be diagonal, and the mass and width parameters in 
Eq. (^) are considered exact values. The values in Eq. (^) lead to the following branching 
ratios, depending strongly on the u decay channel from which one extracts gpujn, 

r,f - - \ _l (0-98 ± 0.21) % g,^^ = 11.7 ± 0.4 GeV-^ ,__. 

b[T ^ujTT u^) - < ^^ g^ ^ Q 23) % g^^ = 15.0 ± 0.1 GeV-i ^"^^^ 



The errors in the branching ratios are dominated by the errors in A and k in Eq. (|24f ). A 
comparison to the measured experimental branching ratios shows that small values for (yfp^^^ 
are excluded, 

i3e^.p. (r- -^ ujTi-Ur) = (1-92 ± 0.08) % (26) 



where we combined the measured branching fractions from CLEO [Q and ALEPH |2^ . 
Thus we will put gp^jn = 15.0 GeV~^ in the following, keeping in mind that the measured r 
decay rate would even require a higher value of gpujn- 

The transition u — > ir'^n' is assumed to proceed through p^u mixing and is written in 
the form 

{7r+{h)7f-{k2)\T\uj{k, A)) = ^-f^BWp{e){h - k2), e^{k, A) (27) 

where gp-^n is related to the decay p° -^ tt+tt^, and the p^u mixing parameter Op^^ is measured 
in e"'"e~ -^ tt+tt^ experiments |]26[| , 

gp^^ = 6.08 , Op^ = (-3.97 ± 0.20) x lO"''^ GeV^ . (28) 

Combining the amplitudes in Eqs. ( PD| , p7D one obtains for the three pion decay mode after 
summation over the polarization A of the intermediate u state, 

(7r-7r-7r+|y"|0) = J2 (vr+7r-|T|cu(j9, X)){uj{p, A) n-\V^\0) ^BW^(s) (29) 

where s = p^ is the momentum transfer and the width in BW^(s) is chosen to be energy 
independent due to its smallness, 

2 
Tfl 

BW^{s) = — ^. — , M^ = 0.782 GeV, T, = 8.4 MeV (30) 

mf^. — s — im,„i ,, 



''LU-^ LU 



With the identity J2x£ij.{p, X)£t{p, X) = —g^" +p^p'' /m? we find the following parametriza- 
tion of the form factor F3, 

(7r-(gi)7r-(g2)7r+(g3)|V'^(0)|0)=ze^°^^gi„g2/3g37^3(si,S2,g') 






X [BW^(si)BWp(si) - BW^{S2)BW p{s2)] . (31) 
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FIG. 3. The r^vrvr mass spectrum from t~ — > Tjir^TT^Ur measured by CLEO [^ (filled circles) 
and by ALEPH |2^ (open circles) normalized to Eg = T{t~ -^ e~Vf,VT)- The solid line shows the 
fit result to Eq. (p3)i the dashed line represents the ?77r7r mass spectrum obtained from e^e~ — > "qin: 
data pl|n]. 



V. VECTOR CURRENT CONTRIBUTION TO r 



-K^Tl^TT Vt 



In the case of the r 



7r°7r°7r v^ 



decay mode we assume that the vector current 



contribution is generated by r^vr mixing in the decay chain r 



rip V-r 



Tl^Tl^Tl V^ 



The r -^ 'qTi^ii u^ decay is allowed to proceed in the Standard Model via a vector 



current induced by the Wess-Zumino anomaly part in the Lagrangian |23] . A normalization 
of the form factor F^ 



(v^T^) jg fixed in the chiral limit and a parametrization of F^^^'^ reads 



{viqi)rf%q2)rf~iqs)\V>^i0m 



«e giQ,g2 /JO'S 7-^3 



i(r77r7r) 



(si,Q' 



^^ m(3m) fQ2\ rp(2ra) ( 



The form and parameters of Tj^™'^(s) are given in Eqs. (|T 



(32) 
For the form factor Ti^"^^ (s) 



the superscript (3m) implies the anomalous transition p —^ rjp —>■ rjirn. In |^,n,5^, a form 
for TJ^^^\s) including p, p' and p" was used, which has been obtained from a fit to (fairly 
poor) e+e~ —>■ tj-r-k data pSJ^ 



However, new measurements for r 



r]7r^7c' 



Ut- have 



become available allowing now for a direct determination of Tj^''")(s) in r decays p5| , |3T|] . A 



10 



e<^T 



direct fit to the differential decay rate for r -^ rin^Ti Vr normalized to Fe = r(r -^ e Ve^- 
as shown in Fig. |^ (solid line) yields for the coefficients ^ and a: 



" ^'^"1 + e + O- 



BWp(s) + eBWp,(s) + crEWp^is) 

^ = -0.22 ±0.03 
(r = -0.10 ±0.01 
xVd.o.f. = 11.0/14 . (33) 

where the masses and widths of the resonances are given in Eq. (^). Again the errors have 
to be considered educated ones, see the remark in Sec. III. The branching fraction that we 
obtain is compatible with the measured decay rate, 



B [t- -^ T]n'^TT-Ur) = (0.14 ± 0.05) % 
Bexp. {r- -> rin\-Ur) = (0.17 ± 0.03) % . (34) 

where we give the weighted average of the experimental branching fractions from CLEO 



1 31 1 and ALEPH |25]. Thus the invariant mass distribution and the decay rate for the 
r^ -^ rjTT^TT"!/^ decay mode are well described by these parameters. On the other hand we 
found that the //tt^tt" invariant mass spectrum in Fig. ^ is only poorly described by the 
TJ^'^^(s) parametrization based on the e"'"e^ -^ 7]nn data (dashed line in Fig. ^. 
For the isospin violating form factor in the three pion decay we deduce the form 



{7f%q,)7T\q2)7T-{q^)\V^{0)\0)=te''''^^q,^q2pq3,F3{suS2,Q 



2\ 



(■»?7r7r)/ ^2\ I?iV^'^) f ^ /n2> 



F,{s^,S2,Q')=e Fr^\s,,Q')-Fr^\s2,Q') , (35) 



where an estimate of the rjir^ mixing parameter e is given by |32 

e = (1.05 ±0.07) X 10^2 _ (gg) 

An additional decay channel would be induced by 7777' mixing with a subsequent r^V 
transition. Experimentally, the decay t~ —>■ rj'iT^TT^Ur has not been observed PBI and thus 



a reliable parametrization of the associated form factor Fg ^'^ is missing. We therefore 
neglect possible contributions from the 77' as an intermediate state. 



VI. NUMERICAL RESULTS 

After having fixed our model for the isospin violating vector current contributions to the 
three pion decay mode, we next discuss numerical effects of this contribution to the decay 
widths and in particular to the structure functions Wb, sign{si—S2)WF, sign{si—S2)WG, Wh 
and Wj. The structure functions Wp and Wg are again combined with an energy ordering 
sign{si — S2) to account for Bose symmetry. 

Let us start with the r^ -^ tt^tt^tt+z/t- decay mode. Fig. ^ shows the resonance structure 
of the pure vector structure function Wb- The (tt^tt^tt^) mass distribution (sohd line) is 
dominated by the two higher radial excitations p'(1465) and p'(1700) of the p resonance 
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FIG. 4. 

Wb in the t~ — 



Invariant mass 
= m{'K^-K^) (dashed) and x 
vr^vr^vr+z^T- decay mode. 
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distributions 



/Q2 



7TT, vr vr vr ' 



(sohd), 



s^ = m['K vr ) (dotted) of the structure function 



in Eq. (pO]). The narrow peak in the y/sl = y^ = m(7r+7r^) distribution (dashed hne) 
shows the dominance of the uj sub-resonance in the vector current. The shape of the y^ 
distribution (dotted hne) is fixed by phase space restrictions and the ^/Q^ and JJl, JJ^ 



distributions through s^ = Q"^ — Si — S2 + 3m^. A comparison of Wb in Fig. ^ with Wa 
in Fig. |I| shows that the contribution to the decay rate from Wb is small compared to the 
axial vector structure function Wa- In fact, using Eq. ([TT| ) we find that Wb contributes 
numerically 0.4% to the decay rate, which is slightly below the branching fraction that has 
been reported by ARGUS [0, namely 0.6%. Due to the large uncertainties in the axial 
vector part, in particular in the ai width, isospin violating effects cannot be seen by a 
rate measurement. One could try to disentangle the structure functions Wa and Wb by 
analyzing the difference in the cos/3 distribution (see [|lO| , pJ] ) (/? denotes the angle between 
the normal of the three pion plane and the direction of the laboratory in the hadronic rest 
frame). However, the sensitivity to the difference in the /? distribution for these two structure 



functions is fairly small |T6| and such an analysis is probably not possible with the current 
statistics. 

Much more promising is an analysis of the vector current contribution through the 
measurement of the interference effects between the vector current contribution with the 
dominating axial vector current contribution, i.e. through a measurement of the structure 
functions sign{si — S2)Wf, sign{si — S2)Wg, Wh and Wj. The three meson and two meson 
invariant mass distributions for these structure functions are shown in Fig. ^ The shape 
of the three pion invariant mass distributions (solid lines) is determined by the interference 
of the tti resonance in the form factors Fi and F2 and the TJ^™^) resonance in Eqs. ( p0l , p2|) . 



The p' and p" peaks are visible in all four structure functions. Similarly, the narrow peaks 
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FIG. 5. Invariant mass distributions x = \fQ^ = m(7r^7r^7r+) (solid), x = a/si^ = rn-(7r^7r^) 
(dashed) and x = ^/s^ = m{7r~TT~) (dotted) of the structure functions sign{si — S2)Wf (a), 
sign{si — S2)Wg (b), Wh (c), Wi (d) in the r" — > i:~ tt~ t:~^ v-j- decay mode. 



around 800 MeV in the m(7r+7r^) invariant mass distribution is a consequence of the inter- 
fering p resonance in Fi and F2 with the product of BWt^(si,2)BWp(si,2) in -^3 as described 
in Eg. (|3TD . The structure functions sign{si — S2)Wf and Wh are the most promising candi- 
dates to extract a vector current contribution in an unambiguously way. An additional scalar 
contribution of the same size as the vector current contribution discussed before would con- 
tribute to two additional structure functions whose angular coefficients are similar to those 
of Wg and Wi |]10[ , and thus a separation of the vector current and such scalar effects might 
be very difficult with presently available statistics in the data (see also |jl6|)- On the other 
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hand, the angular distributions which determine the structure functions sign{si — S2)Wf 
and Wh differ considerably from those originating from possible spin-zero-spin-one interfer- 
ence effects. Any nonvanishing contribution to these structure functions would therefore be 
a clear signal of isospin violation. 

In the decay r^ -^ -k^-k^-k^Vt- we find the effects of isospin violation to be negligibly small. 
Indeed, the contribution to the decay rate from Wb is of the the order 10~^%. Those the 
amplitudes in the invariant mass distributions of Wb are very small when compared to the 
corresponding invariant mass spectra in the decay r~ — » 7r~7r~7r"'"z/^. Even the distributions 
in the interference terms do not have significant amplitudes. We therefore conclude that 
isospin violation in the decay r^ -^ 11^71^71^ Uj. can hardly be measured in presently available 
data. 



To summarize: An isospin violating vector form factor is expected to give a contribution 
of about 0.4% to the decay rate in the r~ -^ tt^tt^tt+z/^ decay mode. Sizable interference 
effects of this vector form factor with the dominating axial vector form factors are discussed 
in detail. These effects could be observed with presently available statistics without recon- 
structing the r rest frame. Any nonvanishing contribution to the corresponding structure 
functions would be a clear signal of isospin violation. The corresponding signals in the 
T^ -^ Tx^Tx^Tx^v^ decay are found to be considerably smaller. 
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